SPECTRAL ASYMPTOTICS FOR THE THIRD ORDER OPERATOR 
WITH PERIODIC COEFFICIENTS 



ANDREY BADANIN AND EVGENY KOROTYAEV 

Abstract. We consider the self- adjoint third order operator with 1-periodic coefficients on 
the real line. The spectrum of the operator is absolutely continuous and covers the real line. 
We determine the high energy asymptotics of the periodic, anti-periodic eigenvalues and of 
the branch points of the Lyapunov function. Furthermore, in the case of small coefhcients we 
show that either whole spectrum has multiplicity one or the spectrum has multiplicity one 
except for a small spectral nonempty interval with multiplicity three. In the last case the 
asymptotics of this small interval is determined. 



1. Introduction and main results 

We consider the self-adjoint operator H acting in L^(]R) and given by 

H = id^ + ipd + idp + q (1.1) 

where the real 1-periodic coefficients p, q belong to the space L^(T), T = M/Z, equipped with 
the norm ||/|| = \f{s)\ds. Without loss of generality we assume that 

[ q{t)dt = 0. (1.2) 
Jo 

A great number of papers is devoted to the inverse spectral theory for the Schrodinger 
operator with periodic potential: Dubrovin [D], Garnett-Trubowitz |GT] . Its-Matveev |IM] . 
Kargaev-Korotyaev |KK] . Marchenko-Ostrovski |M0] . Novikov |No] etc. Note that Korotyaev 
|K3 ] extended the results of |M0] . |GT] . |KK] . for the case —y" + qy to the case of periodic 
distributions, i.e., —y" + q'y on I/^(M), where periodic q G Lf^^(R). 

The results for the operator H are used in the integration of the bad Boussinesq equation, 
given by 

p = ^d^(d^p + 4p^y p = dq, (1.3) 

on the circle, see |McK] and references therein. Here it (or du) means the derivatives of the 
function u with respect to the time (or space) variable. It is equivalent to the Lax equation 
H = HK - KH where K = -d"^ + |p. 

The inverse scattering theory for the self-adjoint third order operator id^ + ipd + idp + q 
with decreasing coefficients was developed in |DTTj . McKean |McKj obtained the numerous 
results in the inverse spectral theory for the non-self- adjoint operator H^, = +pd + dp + q on 
the real line with smooth and sufficiently small p and q. Results for non-self-adjoint operator 
was applied for integration of the good Boussinesq equation. 
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It is known much less about the self-adjoint third order operators with periodic coefficients. 
Even the direct problem is not well developed. 

Results about the spectrum of the higher order differential operators with smooth periodic 
coefficients are given in the book [DSj, see also McGarvey's paper |McG] . The case of the 
even order differential operators with non-smooth coefficients are given in |BK3j . 

The operator H was considered by the authors in [BK4j . In order to describe our main 
results we recall needed results from |BK4j . We consider the differential equation 

iy'" + ipy' + i{py)' + qy = Xy, (t, A) G M x C. (1.4) 

Introduce the 3x3 matrix-valued function M{t, A) by 

V'l V'2 ^3 (^,A), (t,A)GRxC, (1.5) 

y^l+Py^l ^2+P^2 ^l+PVzj 

where Lpi,ip2,'^z are the fundamental solutions of equation (11. 4p satisfying the conditions 

M(0,A) = l3, VAgC. (1.6) 

Henceforth I3 is the 3x3 identity matrix. We define the modified monodromy matrix by 
M(l, A), below it is called shortly the monodromy matrix. The matrix-valued function M(l, ■) 
is entire and its characteristic polynomial D is given by 

D(r,A) = det(M(l,A) -r]l3), (r, A) G (1.7) 

An eigenvalue of M(l, A) is called a multiplier, it is a zero of the algebraic equation D{-, A) = 0. 
If r(A) is a multiplier for some A G C, then r(A)~^ is also a multiplier. Each M(l, A), A G C, 
has exactly 3 (counting with multiplicities) multipliers Tj{X),j = 1,2,3, which satisfy 

Tj{\) = e''''''\l + Oi\z\-')) as |A|^oo. (1.8) 



Henceforth we put 



,^ , , TT 37r 

A3, argAG(--,y 



arg^; G 

The operator H is self-adjoint on the domain 



TT TT 
6' 2 J 



■2n 

w = e 3 . 



Dom{H) =|/ G L'{R) : zif'+pf)' + ipf + qf e L^R), [f' + Pf)' e AocW}- (1-9) 

The spectrum of H is absolutely continuous and satisfies 

a{H) =|a G M : |r,(A)| = 1, some j = 1,2,3| = 61 U ©3 = ©2 = 0, 

where ©^ is the part of the spectrum of H having the multiplicity j = 1, 2, 3. The spectrum of 
H has multiplicity 1 at high energy. Note that the spectrum of the even order operator with 
periodic coefficients has multiplicity 2 at high energy, see jBK3j . Recall that the spectrum 



of the Hill operator —d'^ + g on the real line is absolutely continuous and consists of spectral 
bands separated by gaps. In the case of the "generic" potential all gaps in the spectrum are 
open, see |M0] . |K3j . In the case of third order operators with periodic coefficients there are 
no gaps in the spectrum for all p, q. 

The coefficients of the polynomial D{-,\) are entire functions of A. Due to (11.81) there 
exist exactly three roots ti(-), r2(-) and t^{-), which constitute three distinct branches of some 



SPECTRAL ASYMPTOTICS FOR THE THIRD ORDER OPERATOR 



3 



analytic function r(-) that has only algebraic singularities in C, see, e.g., Ch. 8 in [Foj. Thus 
the function r(-) is analytic on some 3-sheeted Riemann surface TZ. There are only a finite 
number of the algebraic singularities in any bounded domain. In order to describe these points 
we introduce the discriminant p{X), A G C, of the polynomial D{-, A) by 

P = (ri - r^fin - Tsf{r2 - nf. (1.10) 

The function p = p(A) is entire and real on M. Due to McKean |McKj a zero of p is called a 
ramification and it is the branch point of the corresponding Riemann surface TZ. Ramification 
is a geometric term used for 'branching out', in the way that the square root function, for 
complex numbers, can be seen to have two branches differing in sign. We also use it from 
the opposite perspective (branches coming together) as when a covering map degenerates at 
a point of a space, with some collapsing together of the fibers of the mapping. In the case of 
the non-self-adjoint operator H^, the ramifications are invariant with respect to the Boussinesq 
flow and they can consider as the spectral data, see [McKj . 
If p = g = 0, then the function p and its zeros r°'^ have the form 



Po = 64sinh^^^ — sinh^ sinh^ , r"'+ = r^'" = i ^ , n e Z. (1.11) 

2 2 2 \ y/3 J 

We formulate our first results about the zeros r^, n G Z, of the function p (the ramifications). 

Theorem 1.1. i) The function p is entire, real on R, and satisfies 

p(A) = |T(A)|^-8ReT3(A) + 18|T(A)p-27, V A G M, (1.12) 

where 

T(-) = TrM(l,-). (1.13) 

a) Let ©3 be the part of the spectrum of H having the multiplicity 3. Then 

©3 = {A G M : \Tj{X)\ = 1, V J = 1,2,3} = {A G M : p(A) ^ 0}. (1.14) 

Moreover, there exists only a finite number (^ 0) of the bounded spectral bands with the 
spectrum of multiplicity 3. 
Hi) The ramifications as n +oo satisfy: 



^TTTt / 

= rf„ = r°'^ - ^^[Po T \Pn\ + 0{wn) + 0{n-^) 



! p{t)e~'^'''''dt, neZ, 
Jo 



(1.15) 



for some sequence Wn,n G N, such that J2n^i ^^^^ ^ 

Remark, i) Identities fll.l4p show that the endpoints of every spectral interval with the 
spectrum of multiplicity 3 are the ramifications. 

ii) In this paper we analyze the spectrum at high energy using the multipliers. The functions 
Aj = j{Tj + T~^), j = 1, 2, 3, are the branches of the standard Lyapunov function A(A) analytic 
on a 3-sheeted Riemann surface. This function is more convenient for analysis of the spectrum 
at finite energy. Note that identity fll.l4p implies: 

63 = {AgR: A,(A) G [-1,1], Vj = 1,2,3}. 

The graph of a typical Lyapunov function and the spectrum 63 are shown by Fig. [H 
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Figure 1. The function p, the Lyapunov function and the spectrum 63 

Let (A2n)nez be the sequence of eigenvalues of equation (11.41) with the 1-periodic boundary 
condition y{x+l) = y{x),x G M. Let {\2n+i)nGZ be the sequence of eigenvalues of the equation 
fll.4p with the anti-periodic boundary condition y{x + 1) = —y{x),x G M. They are labeling 
(counted with multiplicity) by 

... ^ A_4 ^ A_2 ^ Ao ^ A2 ^ A4 ^ the periodic eigenvalues, 

... ^ A_3 ^ A_i ^ Ai ^ A3 ^ the anti-periodic eigenvalues. (-'-•-'-6) 

If p = g = 0, then the periodic and antiperiodic eigenvalues are given by A° = (vrra)^, n G Z. 

Theorem 1.2. i) The periodic and antiperiodic eigenvalues satisfy 

A„ = (vm)^ — 2pQnn H — ^ as n ^ ±00. (1-17) 

n 

a) The entire function T = TrM(l,-) and the spectrum 63 of the multiplicity three are 
recovered by the periodic spectrum plus one antiperiodic eigenvalue or by the antiperiodic spec- 
trum plus one periodic eigenvalue. 

Remark, i) Asymptotics (ll.lSp of the ramification is sharp, since it is determined in terms 
of Fourier coefficients of p, q. Unfortunately, the asymptotics of the periodic and antiperiodic 
eigenvalues in (I1.17P is not sharp. 
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ii) Knowing the function T and using identities fll.l2p . (12. 6p we recover the functions 
p(A), D{t, A) and all multipliers. 

We consider the operator if^ acting in L^(R) on the domain (11.91) and given by 

H, = id^ + e{ipd + idp + q) (1.18) 

where e G M is a small coupling constant. We have the following result. 

Theorem 1.3. Let p G L^(T) satisfy p{t)dt = 0. Then there exist two functions r^{e), 
real analytic in the disk {\e\ < c} C C for some c> 0, such that r^(0) = and they satisfy 



r+{e)-r-{e)=4:h^e^ + 0{e^) as e ^ 0, (1.19) 

if h>0 
if h<0 



Ar {E),r+{e)) or {r+{e),r (e)) if h > ^g) 



where 



h-lT.(M.-Mi)- (1.21) 



3^V(27rn)2 (2™) 



Remark, i) The functions r^{e) are the (nearest to A = 0) zeros of the function p{X,e). 

ii) Let £ > be small enough. U h > 0, then r^{e) are real and there is a band of the 
spectrum of the multiplicity three. If /i < 0, then r^{e) are non-real and there is no a band 
of the spectrum of the multiplicity three. 

iii) The proof of the theorem is based on the analysis of the monodromy matrix as £ — ?■ 0, 
and identities (I1.12p . (I1.14p . We determine the asymptotics of r^{e) in the form r^{e) = 
r{e) ± 2h2e^ + 0{e'^) as e — > for some function r. This gives the asymptotics (I1.19p . In 
the proof we use the some technique developed for fourth order operator with small p, q, see 
pal[BK2] . 

We describe now the results for vector differential equations and higher order differential 
equations. We begin with the vector case, where more deep results are obtained. The inverse 
problem for vector-valued Sturm-Liouville operators on the unit interval with Dirichlet bound- 
ary conditions, including characterization, was solved by Chelkak-Korotyaev [ CKlj . [CK2j . 
We mention that uniqueness for inverse problems for systems on finite intervals was studied 
in [Ma] . The periodic case is more complicated and a lot of papers are devoted only to the 
direct problem for periodic systems: Carlson [Cal] , |Ca2] . Gelfand-Lidskii [GLj . Gesztesy and 
coauthors [CLj . |CG] . Korotyaev and coauthors |CK] . [BBKj . |K1] . |K2] . etc. The discrete 
periodic systems were studied in [KKul] , |KKu2] . We describe results for the first and second 
order operators with the periodic N x N matrix- valued potential from |CK] . [Klj . |K2] . which 
are important for us. In fact the direct problem is consisted from two steps: 

First step: 

(1) the Lyapunov function on some Riemann surface is constructed and described, 

(2) sharp asymptotics of periodic eigenvalues and ramifications of the Lyapunov function 
are determined, 

(3) multiplicity of the spectrum, endpoints of gaps are the periodic or antiperiodic eigen- 
values or the ramifications of the Lyapunov function are determined. 

The second step is more complicated: 



6 



ANDREY BADANIN AND EVGENY KOROTYAEV 



(4) the conformal mapping with real part given by the integrated density of states and 
imaginary part given by the Lyapunov exponent is constructed and the main properties are 
obtained, 

(5) trace formulas (similar to the case of the Hill operators) are determined, 

(6) global estimates of gap lengths in terms of L^-norm of potentials are obtained. 
Spectral theory for higher order operators with decreasing coefficients is well developed, see 

|BDTj and the references therein. Numerous papers are devoted to the fourth order operators 
on bounded interval: [B], |CPSj . |McLj . [S] etc. The third order operator on the bounded 
interval was considered by Amour |Alj . |A2j . 

Even 4) order operators with periodic coefficients considered in the papers: Badanin- 
Korotyaev jBKl] - pK8] . Papanicolaou [Ell ES], Mikhailets-Molyboga |MM1[[MM2] . Tka- 
chenko |Tk] . see also references therein. The spectral analysis of the higher 3) order 
operators with periodic coefficients is more difficult, than the analysis of the first and second 
order systems with periodic matrix-valued potentials. The main reason is that the monodromy 
matrix for the first and second order systems has asymptotics in terms of cos and sin bounded 
on the real line. The asymptotics of the monodromy matrix for higher order operators has 
additional components in terms of cosh and sinh unbounded on the real line. 

The 2A^-order {N ^ 2) operator with periodic coefficients was considered in [BK3j (the case 
N = 2 see also in |BK1] . |BK2] ) and only the first step was done. The conformal mapping for 
the higher order operator, which is important for the spectral analysis, is not still constructed 
and there are no gap length estimates in terms of the norms of potentials. 

The plan of the paper is as follows. In Sect. 2 we describe the basic properties of the 
monodromy matrix A4. In Sect. 3 we consider the function p and the Riemann surface of the 
multipliers at high energy and prove Theorem 11.11 i).ii). In Sect. 4 we determine asymptotics 
of the ramifications and prove Theorem 11.11 iii). In Sect. 5 we determine asymptotics of the 
periodic and antiperiodic spectrum and prove Theorem 11.21 In Sect. 6 we consider the case 
of the small coefficients and prove Theorem 11.31 The technical proof of the asymptotics of the 
trace of the monodromy matrix Ai is placed in Appendix. 



2. Monodromy matrix 

Consider the unperturbed equation iy'" = \y. It has the solutions e^^'^^ ^* and the multipliers 
have the form e*'^ ^, j = 1, 2, 3, here and below 



u; = e3, z = x + zy = A3 G 5, argA G (^--, yJ, 5 = {2; G C : arg z G (^--, -J }, 

and on the boundary of S we identify each point z = re*2,r > 0, with the point re *6. The 
trace Tr Mq of the unperturbed monodromy matrix Mq is an entire function in A given by 

To = Tr Mo = e'' + e^' + e^^'. (2.2) 

Consider the perturbed equation fll.41) . The matrix- valued function M(t, A), given by fll.Sp . 
satisfies 

M' - P{X)M = Q{t)M, M(0,A) = l3, all (t. A) G M x C, (2.3) 
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where the 3x3 matrices P and Q have the form 

P = \ 1|, Q = \ -p 0|. (2.4) 
\-iX 0/ \iq 

The matrix- valued function M(l, ■) is entire and satisfies 

detM(l,A) = l, M*{1,X)JM{1,X) = J where J= -i \ (2.5) 





for all A G C, see |BK4] . This identity is an odd-dimensional analog of the symplectic property 
of the monodromy matrix for the even order operator [BK3] . Moreover, 

D{t, A) = -T^ + T^T{X) - tT(X) + 1, all (r, A) G (2.6) 

Introduce the simple transformation 

M = ^-^M^, ^-^P^ = izn, Q = '^-^Q^ = — (pQi + -Q2] (2.7) 

6ZZ V Z / 

where = ZU and 



U=^\l u con = {U*)-\ (2.8) 










0^ 




n= 1 





u 









^0 












Qi = I 1 -2u'^ u I , Q2= { tu u u \ . (2.9) 

Using this transformation we rewrite the problem (12. 3 p in the form 

M' -izQM = Q{t,X)M, A^(0,A) = l3. (2.10) 

Identities (12. 7p . (12.91) show that the matrix Q in the right hand side of equation (12.101) 
satisfies Q(t,X) = 0((|p(t)| + |g(t)|)|2;|~"'^) as |A| — )■ 00, while the corresponding coefficient Q 
in (12.31) is only bounded. It is a crucial point for our analysis. 

In [BK4j we proved that the monodromy matrix Ai and its trace T(A) = Tr A^(l, A) satisfy 

\M{1, A) - e'^^l ^ ^e"""^", all |A| ^ 1, (2.11) 

|T(A)| ^ 3e"°+", all A G C, 
|T(A) -To(A)| ^ ^e^o+^ all |A| ^ 1 ^^'^^^ 



where 



y + V^a 



\z\ 



^ = Ibll + Wlh ^0 = Txiax Re{izu^) = Re{izu ) = , — ^ -Zo ^ k|, (2-13) 

j=o,i,2 2 2 

and z = X + iy. Henceforth a matrix A has the norm given by 

\A\ = max{\/E ^ : is an eigenvalue of the matrix A* A}. (2-14) 
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Below we need sharper estimates of the monodromy matrix. Equation (I2.10p and the stan- 
dard arguments provide that the function A^(t, A) satisfies the integral equation 

M{t,X) = e'''^+ [ e''^'-''>^Q{s,X)Mis,X)ds, all (t, A) G R x (C \ {0}), 
Jo 

where the matrix Q is given by fl2.8p . The standard iterations yield 

oo 

M{t, A) = Y,^n{t, A), Mo{t, A) = e^^*^ (2.15) 



where 

Mn{t,\)= f e''^'-''^^Q{s,\)Mn-i{s,\)ds, neN. (2.16) 



Lemma 2.1. The series \2.15\) converges absolutely and uniformly on any compact in C\{0}. 
The matrix-valued function M.{1, ■) is analytic in C\ {0} and satisfies 



|M(1,A)| ^ e^°+", |M(1,A) - ^M„(1,A)| ^ all NeN, |A| ^ 1. (2.17) 

Proof. Identities fl2.16p give 

~ n 

Mn{t,X) = I l[(e"^''+^-''^''Q{tk))Mo{ti,X)dtidt2...dtn, 

0<tl<...<tn<t„ + l=t ^^=1 

the factors are ordering from right to left. Using estimates |e*^*^| ^ e^°* we obtain 

|A^„(t,A)|^— / \Qis)\ds] , all (n,t,A) gNxM+ X (C\{0}). (2.18) 



These estimates show that for each fixed t G M+ the formal series f l2.15p converges absolutely 
and uniformly on any compact in C \ {0}. Each term of this series is an analytic function in 
C \ {0}. Hence the sum is analytic in this domain. Summing the majorants and using the 
estimate \Q{s)\ds ^ ^ we get fl^TTp . ■ 

In the following Lemma (proof in Appendix) we will determine asymptotics of the trace 
T of the monodromy matrix. Introduce the auxiliary function (p which will be used in this 
asymptotics: 

</,(t,A)= Yl a;2(fc+i)giu;^'.gi(a;'»-a;i).i^ (t, A) G M X C. (2.19) 

For p ^ 1, a G R we introduce the real spaces 

P' =[f = UnUl.. ||/||p<Oo|, 11/11^ = ^|/„r<00, 

C ={/ = {fn)n&, 5^(1 + |27rn|2-)|/„|P < ooj. 
We will write a„ = P'^iji) iff the sequence (a„)„ez £ ^a- 
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Lemma 2.2. Let p,q G L^(T). Then the function T satisfies 



T = $o + ^ + ^ (2.20) 



where 

$o = ^e^-^+^, <!>,{X) = -l [\{s,XMs)ds, r^{s) = !\{t)p{t - s)dt, (2.21) 
fc=0 ^ -^0 

anc? t/ie function $ satisfies 

~, ,A 0{\z\-^) as |A|^oo, argAG[-f,^] 

^^^^^''^\e{n) + 0{n-^) as n ^ +oo, \ = [l + 0{n-^)) ' ^^'^^^ 
uniformly in arg A G [— f , 

3. Ramifications 

In this Section we will consider the function p given by fll.lOp . Now we will prove a Counting 
Lemma for the ramifications. Introduce the contours 

Ca{r) = {\e C:\z- a\<r}, a G C, r > 0, 

and the domains 



V+n = \XeC: 



Lemma 3.1. T/ie function p is entire, real on R, and satisfies identity U.1S\) and asymp- 
totics 

p = Po{l + 0{\z\-^)) as|A|^oo, \eC\VJnejJ)n. (3.2) 
ii) For each odd N > large enough the function p has (counting with multiplicities) 2N 

zeros on the domain {|A| < i^)^} and for each \n\ > exactly two zeros in the domain 

T>n- There are no other zeros. 

Hi) The function p has a finite even number 2m ^ of real zeros, counted with multiplicities. 

Proof, i) The standard formula for the discriminant d of the cubic polynomial — + ar^ — 
6r + 1 gives d = {abY — 4(a^ + 6^) + 18a6 — 27 (see, e.g., [Co] . Ch.7.5) which implies that p is 
entire and real on M and satisfies fll.l2p . 
We will show (13.21) . Let |A| — )■ oo. Asymptotics (II. 8p yields 

ri(A)-r3(A) ^ ^ ^ e^^'--')^0{z'^) + 0{z-') ^ 

^iz gia;2z gi(l— tj2)z ' V ' / 

We have 

a/3 

- Rei(l - uj'^)z = Im(l - u^)z = ^(x + V^y) ^ 0. 

Then \e'('^^'^'>\ ^ 1 for all A G C. Moreover, using the standard estimate | sin2;| > ie'^™^' as 
\z — 7Tn\ ^ I for all n & Z (see |PT] . Lemma 2.1), we deduce that 



e 



sm 



2 2 
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as A G C \ Un^i'D_n. Then asymptotics (13. 3p gives 

ri(A)-r3(A) = (e*^-e*-'^)(l + 0(z-i)) as A G C \ U„;,il?_„. (3.4) 
The similar arguments show that 

r2(A)-r3(A) = (e*-^-e*^'^) (1 + 0(^-1)) as A G C \ U„^iP_„. (3.5) 
Furthermore, asymptotics (ll.Sp yields 

ri(A)-r2(A) ^ e'^^-^>0{z-^) + 0{z-^) ^ e^'^-^>0{z-^) + 0{z-^) 



^iz giwz gj(l— Li;)2: ]^ ^ gi(aj— l)z 

The relations 

R n ^ T n ^ J^O,ReA^O 
Re^{l-u:)z = -Ml-u)z = -ix-V3y)i ^^^^^^ 



(3.6) 



give 



e 



^(^-'^^"Kl, allA:ReA<0; |e^(^-^)"| ^ 1, allA:ReA^O. 



Moreover, if A G C \ Un^iVn, then 



(l-co)z 

sm 



> 1 iRc(i(l-a;)2)gilm(l-t^)^ ^ 1 

2 2 



2 

as Re A < 0, and similarly |e*('^~^)^ — 1| = ^ as Re A ^ 0. Then asymptotics (13.61) gives 

ri(A)-r2(A) = (e^^-e*"^) (1 + 0(^-1)) as A G C \ U„5.il)„. (3.7) 

Substituting asymptotics (13.41) . (13. 5p . (13. 7p into (ll.lOp we obtain (13.21) . 

ii) Let ^ 1 be odd and large enough and let A^' > be another odd. Let A belong to 
the contours Co(^),Co(^) and dV^ for all \n\ > Asymptotics ([32]) yields 

|p(A) -po(A)| = |po(A)| ^ - 1 = \poiX)\Oi\zn < |po(A)| 

on all contours. Hence, by Rouche's theorem, p(-) has as many zeros, as po(") in each of the 
bounded domains and the remaining unbounded domain. Since po(-) has exactly one zero of 
multiplicity two at each G Z, and since N' > N can be chosen arbitrarily large, the 

statement i) follows. 

iii) Asymptotics (13. 2p shows that p(A) > on the real line for large |A| > 0. Then p has a 
finite number of real zeros. The function p has even number of zeros in the large disk due to 
the statement ii). The function p is real on R, hence it has both an even number of non-real 
zeros in this disc and an even number of real zeros. ■ 

Recall that the function p has a finite number m of real zeros. Using the results of Lemma 
13.11 we index the zeros r^, G Z, of p by: 
a) labeling of real zeros: 

m even : ^ ^ ... ^ r^m ^ ^ ... ^ ^ r^j", 

m odd : ^ ^ ... ^ ^ ^ ... ^ r^; ^ Tq, 
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b) labeling of non-real zeros: 

m even : < Imrti ^ Imrm , , ^ Imrm , , ^ 

2 2 2 

m odd : < Imr^+i ^ Imr^+i ^ 

2 2 

and 

r-n = rt- (3.8) 

Proof of Theorem 11.11 i). ii). i) The statement is proved in Lemma [3.11 i). 

ii) The first identity in f ll.l4p was proved in |BK4] . We will prove the second one. Let 
Tj = e^''^,j = 1,2,3, where Rekj e (— 7r,7r]. We have ki + k2 + = 0, since T1T2T3 = 1. If 
k i and A G C is not a ramification, then kj{X) 7^ ki{\). Identity f ll.lOp gives 

p = (e'k^ - e*'=2)2(e*'=^ - e*'=^)2(e*'=^ - e'"^^)^ = -64 sin^ sin^ sin^ '^l^, (3.9) 

We have two cases. If A G 63, then the first identity in (11. lip implies that each kj{X) G M, j = 
1, 2, 3, and ([31]) yields p(A) ^0. 

If A G <y{H) \ ©3, then exactly one kj, say ki, is real and k^ = k2 are non-real, since if A G M 
and T2 = e**^^ is a multiplier, then T3 = ^ = e*'^^ is also a multiplier. Thus identity (13. 9p 
implies p(A) = 64| sin p sh^ Im ^2 > 0, which yields the second identity in (II. lip . 

By Lemma [XT] the function p has a finite even number of real zeros. Therefore, there exists 
only a finite number of the bounded spectral bands with the spectrum of multiplicity 3. ■ 

We will show in the following Lemma that the large multipliers in C+ specify by the domi- 
nating multiplier T3 in the sense that these multipliers are zeros of the function 

^(A) = r3(A) - (r3(A))^ (3.10) 

see (I3.15p . Moreover, here we determine the rough high energy asymptotics of the ramifications 
and the multipliers at the large ramifications. 

Lemma 3.2. Let the branches of the multipliers are define by i\1.8\) . Then 

i) The multipliers satisfy 

r-,\\) = ri(A), r,\X) = r.iX) (3.11) 

for all A G A/j and for some R > large enough. 

ii) Let A = and let n — )■ +00. Then 

ri(A) = r2(A) = (-l)"e-^ (l + 0{n-^)), t^{\) = c'tt (1 + 0{n'^)), (3.12) 

ri(A) = r3(A) = (-l)"e^ (l + 0(^-1)) , r2(A) = 6"^ (l + O(n-i)). (3.13) 
Hi) The ramifications satisfy 

= {l ^ 0{n-^)) as n -> +00. (3.14) 

iv) Let A G C+ n for some R> large enough, where the domains A/j = {A G C : |A| > 
R}. Then 

\ is a ramification i^W = 0. (3.15) 
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Proof, i) Recall that r(A) is a multiplier iff r ^(A) is a multiplier. Asymptotics (11 .Sp and 
identity uj = uj^ give 

rri(A) = e-(l + Od^r^)), r^- 1(A) = e-'^(l + 0{\z\-^)), r,\X) = e-^(l + 0{\zn) 

as |A| — i- oo, which yields (13. lip , 
ii), iii) If A = r^, then Tj{X) = rfc(A) for some j, k = 1,2, 3. Let n — )■ +oo. By Lemma I3.H 



/ -i-xi .•7r/27m ^, /27m 



e^6 (^+5„) where < Asymptotics (US]) gives Tj{X) = 0(e~^), j = 1,2, 
and 

r3(A) = e'7f+^"(l + 0(n-i)). (3.16) 

These asymptotics show that T3(A) 7^ Tj(A) for j = 1,2 and for all n G N large enough. Then 
ri(A) = r2(A), which implies the first identity in (I3.12p . 
We will prove (13.140 . Let A = and let n — > +00. Then 

/2iTn ^. /27m ^, 

(_ + o(„ '))=.(_ + o(„ 

and asymptotics (II. 8p yields 

ri(A) = e"(w+'^")(l + 0(n-i)), r2(A) = e"'(^+'") (l + 0(^-1)) . (3.17) 

Substituting these asymptotics into the first identity in (I3.12p and using cu — cu'^ = we 
obtain e'^" = l + 0{n~^). Then 6n = 0{n-^), which yields fIXTip . 

Substituting (5„ = O^n"^) into (13.161) . (I3.17P we obtain asymptotics (I3.12j) . These asymp- 
totics together with fl3llD yield fl3J3|) . 

iv) Let A G C+ n Ar be a ramification. Identities (13. lip and (I3.12p yield f^^{\) = T2(A) = 
ri(A). Using the identity 

ri(A)r2(A)r3(A) = 1 (3.18) 
we obtain ip{X) = 0. Conversely, let ip{X) = for some A G C+ fl A/j. Then identity (13. lip 
implies T3(A) = T2~^(A). Identity (I3.18P yields ti(A) = T2(A), therefore A is a ramification. ■ 

Remark. 1) We describe the surface of the multipliers for the case p = q = 0. We have 

1 

r°(A) = e*^^ in this case, then our surface is the 3-sheeted Riemann surface of the function 
As. We have two parametrization of this surface. 

First parametrization. We construct the standard (escalator type) parametrization TZ^ of 
the surface of the function As if we take 3 replicas TZpj = 1,2,3, of the cut plane C \ iM- 
and join the edge of the cut on the sheet TZi with the edge of the cut on the sheet 7I2, the 
edge of the cut on the sheet 7^2 with the edge of the cut on the sheet TZs, and the edge of the 
cut on the sheet 71^ with the edge of the cut on the sheet TZi, in the usual (crosswise) way, 
see Fig. [2] a). 

Describe the function t°(A). Consider the sectors Sj = u^^^S,j = 1,2,3, on the 2;-plane, 
where S is given by (12. ip . In each of the sectors the function A^, and then t°(A), is univalent. 
Moreover, for each j = 1,2,3, the function r°(A) satisfies: t°(A) = t°(A) = e*"^^ ^ in the 
sector ijj^~^S. The function X = maps each of the sectors Sj,j = 1, 2, 3, onto the sheet 7^° 
of the surface TZ^. For each j = 1, 2, 3, the function t'^(A) satisfies: t°(A) = Tj'(A) on the sheet 

Second parametrization. In order to consider below the perturbed case it will be convenient 
to use the other parametrization TZ^ of the Riemann surface of the function As, see Fig. [2]b). 
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Figure 2. The Riemann surface of the function A3 , which coincides with the Riemann sur- 
face of the muhiphers for the case p ~ q = 0, in a) the standard (escalator type) parametriza- 
tion, b) the parametrization convenient for the considering of the perturbed case 



We take 3 replicas of the cut plane 7^? = C \ «+, = C \ and 7^^| = C \ We obtain 
the Riemann surface TZ^ by joining the edges of the cut iM+ on TZ^ and on the edges 

of the cut z]R_ on 7^° on TZ'^ in the usual (crosswise) way. 

If we deform continuously the surface TZ^ so that the right half-plane Re A > of the sheet 
7Zl and the right half-plane of the sheet 7^° swap places, then we obtain the surface 7Z^. The 
function t^{X) satisfies: 

a) r°(A) = r3°(A) on the sheet TZf, 

b) T°(A) = T°(A) on the left half-plane ReA < of the sheet TZ^ and on the right half-plane 
ReA > of the sheet 7^^; 

c) T°(A) = T2 (A) on the left half-plane of the sheet 7^2 ^^e right half-plane of the sheet 

7^?. 

2) Consider the perturbed case. The surface 7Z of the multipliers for large |A| is close to the 
surface TZ , see Fig. O Let = (^)3, where N is given by Lemma [3.11 11). Then G Vn for 

each |?7,| > ^Y^- Describe the surface TZ for |A| > R. We take 3 replicas TZi,TZ2 and TZ^ of the 
cut plane: 

7^l n A/j = Ah \ u„6Nr„, 7^2 n a^ = a^ \ u„6zr„, 7^3 n a^ = A^^ \ u„,eNr_„, 

where = {A G C : |A| > R}, r„ = [r^.i,?"^] C Au,n G Z, are straight lines. The surface 7Z 
at large |A| is obtained by joining the edges of the cuts r„ on TZi with the edges of the same 
cuts on 7Z2, the edges of the cuts r_„ on 7Z2 with the edges of the same cuts on TZs in the 
crosswise way. 
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Figure 3. The Riemann surface TZ of the multipHers and the ramifications for the case of 
the coefficients po + p, q, with small p, q and constant po- 



4. ASYMPTOTICS OF THE RAMIFICATIONS 

Introduce the entire function ^ by 

e(A) = 4r(A)-T'(A), AgC. 
Lemma 4.1. Let A = r^. Then the function ip{X) = ts{X) — t|(A) satisfies 

V^(A) = Q + 0(n-i))(e(A) + 0(e~^)) as n ^ +oo. 

Proof. Identity (12. 6 p provides 

D(r(A), A) = -r^(A)(r(A) - T(A) + t-\X)T(X) + t~\X)), 

DiT(X),X) = -r{X){T'{X) -t{X)T(X) + T{X) -r-\X)) 
for all AgC. Let A = and let n +oo. Identity (14. 3 p and asymptotics fl3.12p imply 

D(r3(A),A) = -r|(A)(r3(A) -T(A) + 0(e-^)). 
The identity D{t3{X),X) = yields 

r3(A) = r(A) + 0(e-^). (4.5) 



(4.1) 

(4.2) 

(4.3) 
(4.4) 
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Identities flXT^ . flO|) and asymptotics flXT^ . fH3]) give 

D{n{X),\) = -Ts(\){ni\) -Ts(\)T(\) + T{\) + 0{e-^^)) 

= -r3(A)(2T(A) - Ts(\)T(\) + 0(6"^)). (4.6) 
Asymptotics fl3.12p . fl3.13p and identity T = ti + T2 + give 

r(A) = e^(l + 0(n-i)), T(A) = 2(-l)"e^(l + 0(n-i)). (4.7) 
Asymptotics (gT]) and the identity U(r3(A), A) = yield 



r^(A) = ^2;^ + 0(e-^). (4.8) 
T (A) 



Asymptotics (143|) . (gl8|) give 



m = T,{X) - rliX) = T(A) - + 0{e-Ts). 

T (A) 

Then asymptotics (gZD yields fO]) . ■ 
Below we will often use the following simple result. 

Lemma 4.2. Lei / G C([0, 1]) and leta>0,(3 e R. Then 

^i-a+m.uf^^^^^^ fiO) + o{l) 3.^+00. (4.9) 

[a — ip)x 

Proof. We have 

ei-^+iP)^^f(y\du = ^/^^ ^ + / e(-°+^^)^"fc(u)ciM, all x > 0, (4.10) 

(a - z/3)x Jo 

where k{u) = f{u) — /(O). Furthermore, 







A 



max[o,5] |/i;(m)| + e~°'"' max[o,i] \k{u)\ 



^max\k(u)\ f e-"^"dM + max|A;(M)| /" e-"^"(iM ^ 
[0,5] ' "io [0-1] A 



ax 



for any 5 G (0,1). Let x — )■ +cxd and let 6 = Then A = o{x~^). Substituting this 

estimate into f l4.10p we obtain fl4.9p . ■ 

The following Lemma gives the asymptotics of the function ^, given by (14. ip . at the large 
ramifications in C+. 

Lemma 4.3. Let p,q E L^(T). Let X = = z^,( = i'^'^z, and let n — )■ +00. Then the 
function given by ^.1]), satisfies 



e(A) = 4e' 



^ 3C 



(4.11) 
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Proof. Asymptotics fl2.20p yields 

e = eo + | + | (4.12) 

where 

^o(A) = 4$o(A) - $o(A), 6(A) = 4$i(A) - 2$o(A)$i(A), (4.13) 

e(A) = 4$(A) - 2$o(A)$(A) + 0(e""n-^). (4.14) 
Asymptotics fl3.14p yields ( = ^ + 0{n^^). Substituting this asymptotics into identity 
(KTDf we get 

$o(A) = e^-t' + 0(e^^), $o(A) = 2e^-t cos(^ + ^) + 0(e~^). (4.15) 

Substituting fl2:22|) . fl4J5|) into (Oil) we obtain 

^(A) = {i\n) + 0{n'^)) (4.16) 

Asymptotics fl4.15p give 

{„(A) = ie<-9 .in' + J^) + 0(e-?l). (4.17) 

Assume that 

6(A) = ^u;V(|p„p + £l(n) + 0(n-2)). (4.18) 
Substituting asymptotics flirT7|) . f l4?T8|) into f H?T2|) we obtain 

which yields (14. lip . 
We will prove fliTTSj) . Substituting identity fl2:2T]) into fliTTaj) we obtain 

2 /"^ 

^^(A) = -- / a(M,A)r7(n)ci?x (4.19) 



_4n 

72 



9 Jo 

where ^ 

r]{u)= / p{t)p{t - u)dt, a(M, A) = 20(m, A) - $o(A)0(m, A), (4.20) 

is given by f l2.19p . Identities f l2.19p imply 

0(t, A) = e^e("'-^)^* + we^e^"-^)^* + w^g-'Ce^"""')'^*, all (t, A) G M x C. 
Using uj — io'^ = za/3 we have e'^^'^e*^'^"'^^)'^* = 0(e~ vl) and then 

0(t, A) = (e('"'-i)'^* + we^'^-i)^* + 0(6"^™)) (4.21) 
uniformly on t G [0, 1]. Moreover, 

A) = e-^'^e("'-^)^* + u;e-"^e("-^)'^('-*) + w^e-^'^e^"'"")'^*, all (t, A) G M x C. (4.22) 
Asymptotics fl4.15p implies 

$o(A) = 2(-l)"e2 (1 + 0(^-1)). (4.23) 
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Relations ( K^ . yield 

$o(A)0(t, A) = 2(-l)V(e^^^e("'-^)^* + c^e-^^^e^-^^^^^-*) + u'e'^^^e^^'-^^^') (l + 0{n-')) 

= 2e^ (e('^'-i)^* + a;e(--^)^(i-*) + u'e^^"-'^^<') (l + 0{n~')) (4.24) 

uniformly on t G [0, 1]. Substituting identities (SSI]), into we obtain 

a(t,A) =2e^((c^e('^-^)^*-cue('^-^)^(i-*)-(^V('^-"')^*)(l + 0(n-i))+e(^^ (4.25) 
uniformly on t G [0, 1]. Using the identity ri{l — t) = rjit) we obtain 

g(— i)Ct^(^)^i = I e^-'~m^-t)r^^t)dt. (4.26) 



Substituting f l4.25p into f l4.19p and using f l4.26p we obtain 

^,(A) = ^u^e< (^e("'^")^*(l + O(n-i)) + e("-^)^*0(n-^) + e^^'-'^'^'0{n-')^r]{t)dt 

+ 0(n-^)) + e-(^-')™*0(n-i) + e-(^+*)""*0(n-^))r/(t)cit. 



= -u; / I e 

Asymptotics (gl]) and identity e-''^''"-^7]{t)dt = |p„|2 give fHl^ . ■ 

Now, using the results of the previous Lemmas, we will determine the asymptotics of the 
ramifications. 

Proof of Theorem 11.11 iii). Identities = are proved in (13. 8p . Let A = and let 
n +00. Then (1334]) implies C = i^^'^z = ^ + where 5„ = 0{n-^). Relations (1315|) . 

give ^(A) = 0(e ^). Substituting (14.111) into this asymptotics we get 



sin^ 



Recall the estimate Kw"^ + 5^)2 — ty| ^ |£| for all w,e e C (see, e.g., jCKPj . Ch.4.5). Substi- 
tuting 

w=\pnl e= U\{n) + 0{n-^)f' = e,{n) + 0{n-^), 
into the last estimate we obtain 



'\Pn? + ^\ {n) + 0{n-^)) " = \pn\+ i\ (n) + 0{n^). 
^ 2 / 2 

Asymptotics (14.271) gives 

vst™ 3™ V 2 / 

which yields asymptotics (I1.15p . Theorem II. II is proved. ■ 
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5. The periodic spectrum 



In this Section we consider the periodic and antiperiodic eigenvalues labehng by fll.lGp . 
These eigenvalues are zeros of the entire functions D(±l, •), where D is given by (11.71) . Iden- 
tities (12. 6|) give 

D(1,A) = 2iImr(A), D(-1,A) = 2 + 2Rer(A), all A e R. (5.1) 
If p = g = 0, then the functions D(±l, A) have the form 



z zu zu 
Dq{1, A) = -8i sm - sm — sm — 



, . z zu zu 

L'o(— 1, A) = 8 COS - COS — COS — (5.2) 



Now we will prove a Counting Lemma for the periodic and antiperiodic eigenvalues. Introduce 
the domains 



ICr 



|a G C : 1^ - 7rn| < /C_„, = E C : \uz + nn\ < n eN. 



Lemma 5.1. i) The functions D{±1,X) as |A| — t- oo satisfy 

D{1, A) = Do(l, A) (1 + 0{\z\'')) , A G C \ U„ez/C2n, 

D(-l,A) = Do(-l,A)(l + 



A G C \ U„ez/C: 



2n+l- 



(5.3) 
(5.4) 



a) For each odd N > uq for some no ^ 1 the function D{1,-) has exactly N zeros, counted 
with multiplicity, in the disk {A : |A| < (vrA^)^} and for each even n : \n\ > N exactly one 
simple zero in the domain /C„. There are no other zeros. 

Hi) For each even N > no for some uq ^ 1 the function D(— 1, ■) has exactly N zeros, 
counted with multiplicity, in the disk {A : |A| < (ttA^)'^} and for each odd n : \n\ > N, exactly 
one simple zero in the domain K,n- There are no other zeros. 



Proof. We consider only the function -D(l, ■). The proof for D(— 1, ■) is similar, 
i) Identities (15.11) and estimates (I2.12p imply 

|D(1, A) - 1^0(1, A)| ^ 2|T(A) - To(A)| ^ ?^e^«+^ 



all lAI > 1. 



(5.5) 



Substituting the estimate | sinw| > 
the relations 



1 „| Imw)| 



as \w — Tm\ ^ ^ for all n G Z into (15. 2 p and using 



Im2;| + I Im^wl + I Imzw 



\y\ + 



y = 2^0 



we obtain 



|Do(l,A)| =8 



sm 



zu 



sm 



zu 



sm ■ 



> -62 



|(| Im^l+I Im zuj\+\ Im zlu^\) 



(5.6) 



for all A G C \ U„ez/C2„. Estimates ([S3]) and ([5SD yield (K^ . 

ii), iii) Let ^ 1 be odd and large enough and let A^' > be another odd. Let A belong 
to the contours Co(7rA^), Co(vrA^'), C7m(|), \n\ > N, n is even. Asymptotics (15. 3p yield 



|D(l,A)-Do(l,A)| = |Do(l,A)| 



Do{l,X)0{l) 



< \Doil,X)\ 



on all contours. Hence, by Rouche's theorem, D{1, ■) has as many zeros, as Do{l, ■) in each of 
the bounded domains and the remaining unbounded domain. Since Dq{1,-) has exactly one 



SPECTRAL ASYMPTOTICS FOR THE THIRD ORDER OPERATOR 



19 



simple zero at each = (vm)^, n e Z, and since N' > N can be chosen arbitrarily large, the 
statement for D{1,X) follows. ■ 

The following Lemma shows that the periodic and antiperiodic eigenvalues A„ at high energy 
are zeros of the function — 1. The other multipliers remote from the unit circle at large real 
A. Nevertheless, these eigenvalues can specify by the multiplier T3 as well as the ramifications. 



Moreover, we will show that the function — (—1)" is exponentially close to zero at the 



T(A) 

large eigenvalue A„, and determine the rough high energy asymptotics of these eigenvalues. 
Lemma 5.2. Let n — )■ +00. Then the periodic and antiperiodic eigenvalues satisfy 



:-ir = -i(An) = ?7t4 = 2^(1 + 0(e-^-)), (5.7) 
A„ = (7^)3(1 + 0(n-2)). (5.8) 



Proof. Let A = Xn,n — ^ +00. Lemma [5.11 gives z = vm + 5„ G M and < |. Asymptotics 
dM]) gives 

r3(A) = 0(e^-^). (5.9) 
Then |t3(A)| 7^ 1 and identity (13.111) implies T2(A) = rg ^(A). Then 

^"'^ -"^^^^-^^(aRa)-^' 

which yields the first identities in (15. 7p . 
Estimates f l232|) yield 

T(A) = ro(A)(l + 0{z-^)) = {e'' + e'"" + e'"'")(l + 0{z-^)) = 0(e^™). (5.10) 
Substituting (15. 9 p and (I5.10p into (12. 6p we obtain 

D{r,{X), A) = r|(A) {-r,{X) + T(A) + 0(1)). 
The identity D{t3{X), A) = and (KWf imply 

r3(A) = T(A) + 0(1) = T(A)(1 + O(e-^n), 

which gives the last asymptotics in (15. 7p . 

Asymptotics (II. 8p implies (—1)" = ti(A) = e"(l + 0(n^^)). Substituting z = irn + 6n into 
this asymptotics we obtain e*^" = 1 + 0{n^^). Then 5„ = 0(n^^) and z = '7in + 0{n^^) which 
yields (ES}- ■ 

Proof of Theorem 11.21 i). Let A = A„,n +00. Asymptotics (15. 8p shows that z = 
vrn + (5„ e M, Sn = 0(n-^). Asymptotics (12T2UD yields 

r(A) = e-""-a (1 + 0(e-#™)) + + £(;p) . (5.11) 

Identity fl2J9|) gives 

0(t,A) = e^^^'(e^('^-^')"*+we^(^-^')"*+u;2e*(^-"')(^+*)") = e^"'^'(e-^"*+we-^^~'^"*+0(e-^"")) 
uniformly on t G [0, 1]. Substituting this asymptotics into (I2.2ip and using (14. 9 p we obtain 

$i(A) = ( J (e-^^" + a;e-^^"'*^")r/(M)rfM + 0(6"^"")) = e*^"'o(n-i). 
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Substituting the last asymptotics into (15. lip we get 

T(A) = e'^"'+l^§^ (1 + o(n-3)) , T(A) = e"'^""^ (l + o(n-3)) . 
Substituting these asymptotics into (15. 7p we obtain e~*^~ = (—1)" + o{n~^), which gives 



Then 



g 3(^n + 5„) = 1 + 0(n"n 



e 2po . _3s 2po 



which imphes (I1.17P for n — )■ +oo. The asymptotics for n — t- — oo can be obtained by 
substituting —t instead of t in equation (11.41) . ■ 
We need the following Hadamard factorizations of the functions D(±l, ■). 

Lemma 5.3. The functions D{±1,X) satisfy 

D(l,A)=z(Ao-A) n (5-12) 

rieZ\{0} 2n 

^(-l^A) = 8n%^, (5.13) 

uniformly on any bounded subset ofC, where A° = {7m)^,n G Z. 

Proof. We consider the functions D{1,-). The proof for D{—1,-) is similar. Asymptotics 
(I1.17P show that the infinite product in (15.121) converges uniformly on any bounded subset of 
C to the entire function of A, whose zeros are precisely X2n,'n G Z- Identity (15. 2 p yields 

Do(l,A) = -^An%^- (5.14) 

Asymptotics (15. 3 p and identity (15. 2p show that the entire function D{1, X) has order 1. More- 
over, its zeros have asymptotics (I1.17p . Then 

D{1, A) = ^e^^+^(Ao - A) n ^'xo~ ^ ^ for some A,BeC. (5.15) 

Identities flHAijl . fl57[5|l yield 

log = AA + 5 + log ^ + log (5.16) 

Let A — )■ +ioo. Asymptotics (15. 3p implies log ^^^^j^ = 0(|2;|~^). Moreover, we have 



o(l). 



Identity f l536|) gives A = 5 = 0. Identity f l515|) gives fl5:T2|) . ■ 

Now we will prove the results about the recovering of the function p and the spectrum cr{H) 
by the periodic (or antiperiodic) spectrum. 

Proof of Theorem 11.21 ii) . Let {A„, n even} be the periodic spectrum and A,,, be the antiperi- 
odic eigenvalue. Identities (I5.12p . (15.11) give the function D{1,-) and ImT(A) = — |D(1,A). 
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Then we reconstruct ReT(A) up to some constant. This constant, and then the function 
T, can be determined from the identity ReT(A*) = 1,A*) — 1 = —1. Identity fll.l2p 

provides the discriminant p. Identity fll.l4p gives the spectrum ©3 of the multiplicity three. 

Using the similar calculations we recover the function T and the spectrum of the multiplicity 
three by the antiperiodic spectrum and one periodic eigenvalue. ■ 

6. Small coefficients 

We prove Theorem II. 31 Here we use some methods developed for fourth order operators 
with the small 1-periodic coefficients [BKl] . |BK2] . We consider the equation 

iy'" + e{ipy' + i{pyy + qy) = Xy, A G C. (6.1) 

In this case the matrix-valued function M(t, A) = M{t, X,e), (t, A,^) G M x C x M, given by 
(11. 5p is a solution of equation 

M' - P{X)M = eQ{t)M, M(0,A,e) = l3 (6.2) 

where the 3x3 matrices P and Q are given by (12. 4p . 

Consider the case e = 0. The matrix-valued function Mo(t, A) = M(t, A, 0) has the form 
Mq = e*^. Each function Mo{t,-),t G M, is entire. Eigenvalues of the matrix Mq have the 
form e*^*, e*^^*, e*^^^*, since eigenvalues of the matrix P are given by iz, iuz, iu'^z. Estimates 

\Mo{t, A) I ^ e'«l*l, all (t, A) G M x C, (6.3) 

where zq is given by (12.131) and a matrix norm is given by (12.141) . 

Consider the case e ^ 0. The solution M(t,X,e) of problem (16. 2p satisfies the integral 
equation 

M{t,X,e) = Mo{t,X) +e [ Mo{t - s, X)Q{s)M{s, X,e)ds. (6.4) 

Jo 

The standard iterations in (16.41) lead to the standard series 

M(t,A,£) = ^£"M„(t,A), Mn{t,X)= [ Mo{t-s,X)Q{s)Mn-i{s,X)ds, n^l. (6.5) 
n^o Jo 

We need the uniform estimates of the monodromy matrix M(l, A,£:). 

Lemma 6.1. For each t G M the series Ii6. 5\) converges absolutely and uniformly on any 
bounded subset o/C^. Each matrix-valued function M{t, ■),t G [0, 1] is entire in {X,e) G 
and satisfies: 

N-l 

|M(l,A,e)| ^ e^«+", |M(1, A, e) - ^ M„(l, A, e)| ^ |e|^x^e"o+l"l" (6.6) 

n=0 

for all {N, A, £:) G N X where x = + . 

Proof. Consider the case t ^ 0. The proof for t < is similar. Identity (16.50 gives 

/n 
Yl (Mo(tfc+i - tk, A)Q(4))Mo(ti, X)dtidt2...dtn, (6.7) 

0<tl<...<tn<t„ + l=t 



22 



ANDREY BADANIN AND EVGENY KOROTYAEV 



the factors are ordering from right to left. Substituting estimates f l6.3l) into identities l p.7\i we 
obtain 

pZot , i-t n 

\Mn{t,\)\ ^ —i^J \Q{s)\dsj , all (n,t,A) G N X M+ X C. (6.8) 



These estimates show that for each fixed t ^ the formal series (16. 5 p converges absolutely 
and uniformly on any bounded subset of C^. Each term of this series is an entire function of 
{X,e). Hence the sum is an entire function. Summing the majorants and using the estimate 
Jo \Qis)\ds ^ X we obtain (16. 6p . ■ 

We introduce the discriminant p(A, e), (A, e) G C^, of the polynomial det(M(A, e) — rl^) by 
identity f lTTU]) . 

Lemma 6.2. i) The function p(A,e) is entire in and satisfies: 

p{X,e) = po{X){l + 0{e)) as e ^ 0, (6.9) 

uniformly in X on any bounded subset ofV = C\ Ung^I^n, are given by liS. 

a) Let \e\ < c for some c > small enough. Then the function p{-,e) has exactly two 
zeros, counted with multiplicities, in each domain T>n,n G Z. There are no other zeros. In 
particular, the function p{-,e) has no any real zeros in the domain \X\ ^ 1. 

Proof, i) Identity (I1.12p and Lemma 16.11 show that the function p(A, e) is entire. Estimates 
(16. 6 p give M{l,X,e) = Mo(l,A) + 0{e) as e — )■ uniformly in A on any compact in C. 
Let A G P. Then all eigenvalues e*^, e*"^^, e*'^ ^ of the matrix Mo(l,A) are simple and the 
standard matrix perturbation theory (see, e.g., |HJj . Corollary 6.3.4) gives that the eigenvalues 
Tj{X,e),j = 1,2,3, of the matrix M(A,e) satisfy 

Tj{X,e) = e''''~'' + 0{e) = e''^'~''{l + 0{e)) as e^O 

uniformly in A on any bounded subset of V. Substituting these asymptotics into (ll.lOp and 
using the identity 

we obtain (16. 9p . 

ii) Using asymptotics (16. 9p and repeating the arguments from the proof of Lemma 13.11 ii) 
we obtain the statement. ■ 
Introduce the entire functions T„(A) = TrM„(l, A),n ^ 0. Estimates (16. 6p imply 

T(A,e) = TrM(l,A,e) = To(A) +eTi(A) +e2r2(A) +e3T3(A) + 0(£^) as e^O (6.10) 

uniformly in A on any compact in C Below we will use the following relations. 

Lemma 6.3. Let p G L^iT) satisfy p{t)dt = 0. Then the functions Ti and T2 satisfy 

Ti = 0, (6.11) 

ReT2(A) = -3/i(l + 0(A)) as A^O (6.12) 

where h is given by M.21\) . 
Proof. Identity (16. 5p implies 

Ti = Tr Mi(l, ■) = Tr f e^^''^^ Q{t)e'^ ds = Tre^ f Q{t)dt = 0, 

Jo Jo 
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which yields fl6.1ip . Moreover, 

T2 = TrM2(l,-) = Tr /" dtf e^^-'+'^^Q{t)e^'-'^^Q{s)ds. (6.13) 

Jq Jo 

We rewrite identities (12.41) for the matrices P, Q in the form 





1 






^0 





























0) 




u 








Using the identities 



1 t 

1 



1 t 



.tPo -Ti^ + tPo + -P^ = I 1 t 



2 

we obtain 

e*^Wg(,) = e*^°g(s)(l3 + 0(A)) = {-pis)Kr{t) + zg(s)K2(t))(l3 + 0(A)) 



as A — 7- uniformly on {t, s) G [0, 1]^ where 



Ki = e*^°Po* = I 1 t , i^s = e*^oPi = t 
.0 10/ \1 

Using the identities 

Tr Ki(l - u)K^{u) = 2(1 - u)u + ^ ^ ^ + y = - + m(1 - m) 

Trir2(l-n)ir2(^) = ^^^|£^, 

we obtain 

ReTre(i-*+^)^Q(t)e(*-^)-f'g(s) = J(t, s)(l + 0(A)) 
as A — 7- uniformly on (t, s) G [0, 1]^ where 



J(t, s) = p{t)p{s) + n(l - n)) - g(t)g(s)^^^^, n = t - s. 
Substituting this asymptotics into (16.131) we obtain 

ReT2(A)= f dt f J{t,s)ds{l + 0{\)) = ]- [ dt [ J{t, s)ds{l + 0{\)) 
Jo Jo 2 Jq Ji_i 

as A — )■ 0, since J(t, s) = J{s,t — 1). Using the simple identity 

fdt f g{t-s)f{t)f{s)ds = Y,\h\% 

for all / G Li(0, 1),^,^' G L2(0, 1),^(0) = ^(1) = 0, and identities 

f m(1 - M)e-'2™"rfu = /" ^2(1 - u)='e-'2™"dM = 

Jo 2(7^)2' ^ ^ 2(7^)4 

for all n 7^ 0, we obtain (I6.12p . ■ 



24 



ANDREY BADANIN AND EVGENY KOROTYAEV 



Now we will determine the asymptotics of the first spectral interval of multiplicity 3 for the 
small coefficients. 

Proof of Theorem 11.31 Identities (12 ■2p imply 

To(A) = e*^ + e'"^ + e'"'^ = 3 - ^ - ^ + 0(A3) as A ^ 0. (6.14) 

Recall that the functions T{X,e), p{\,e) are entire in {X,e) G C^. Let the entire functions 
a{X,e),b{X,e) and the numbers bj,j = 2,3, be given by 

T{X,e) = 3 + a{X,e) +ib{X,e), b{X,e) = lmT{X,e), all (A,£)GM^ 6^ = ImT,(0). 

Substituting relations flOTj) . f l02|) . IKl^ into f lOOj) we obtain 

T(A, e) = 3 - ^ - ^ - (3/1 - zb2)e' + e'Ts{0) + 0{X') + 0{Xe') + Oie') (6.15) 

as (A,e) — !■ (0,0). Asymptotics (16.151) gives 

a{X,e) = -3h6^ + 0{X^) + 0{X6^) + 0{6^), b{X,e) = fi{X,6) + 0{X^) + 0{Xe^) + 0{e'^) (6.16) 
as {X,e) — )■ (0,0) where 

/^ = -^ + ^, r{E)=2b2E'' + 2bsE^ (6.17) 

Identity (I1.12p gives 

p = a3(a + 4) + 6^ (l08 + 2(a + 18)a + 6^) . (6.18) 
Substituting asymptotics (16.161) into (I6.18P we obtain 

p(A, e) = /(/i, e) = 108 (^fi^ - /iV + 0(/i^) + 0(/iV) + 0(/i£^) + 0(£^)) 

as (/i, £:) — 7- (0,0) where 

A = r-2/i (6.19) 
and f{ii,e) is entire in {fi,e). Introduce the new variable u hj fi = ue^. Then 

/(Me^e) = 108e^E(M,e), E(m, e) = _ /^s + c»(£) (6.20) 

as £ — )• uniformly on any compact in C, the function E{u,e) is entire in {u,e). 

Consider the equation E{u,e) = 0. Let h ^ 0. Using ^_E(±/i2,0) 7^ and the Implicit 
Function Theorem we deduce that there exist two real analytic functions u^{e) in the disk 
{\e\ < c} for some c > such that each u'^i^e) is a zero of the function E{-,e). Asymptotics 
IKW) yields 

u^{e) = ±0 +0{e) as £-^0. 

Substituting /i = u'^{e)e^ into the identity (I6.19P we deduce that there exist two real analytic 
functions r^(e) in the disk {|£:| < c} such that 

r^{e) = r{e)±2hle'^ + 0{e^) as e -> 

which yields (I1.19p . 

Consider h > 0. Let £ > be small enough. Then r~{e) < r{e) < r^{e). Asymptotics (I6.19P 
shows that p{r{e),e) = f{0,e) < 0. Then p{-,e) < on the whole interval {r~ (e) , (e)) and, 
by Lemma [6.21 iii). p{-,e) ^ out of this interval. Then, due to (I1.14p . the spectrum of in 
this interval has multiplicity 3 and the other spectrum has multiplicity 1. The proof for the 
case e < is similar. 
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If /i < 0, then, by Lemma [6.21 iii). the function p has no any real zeros and p{-,e) > on 
the whole real axis. Hence all the spectrum has multiplicity 1. ■ 

7. Appendix 

In this Section we will prove Lemma [2.21 Introduce the functions ejk{t, A) by 

ejkit, A) = e'^'~''e'^^''~'-'^'~'^'\ all {t, A) G M x C, j,k = 1, 2, 3. (7.1) 
Lemma 7.1. Let f,g& L^(T). Then all functions ajk, 1 ^ j < k ^ 3, given by 



ajk{X) = / ejk{t,X) / f{u)g{u-t)dudt, (7.2) 

satisfy asymptotics 

.0/ ^d^l"') as\X\^oo, argAehf,^] 
"^■'^ ^ \e{n) + 0{n-^) as +00, X = -z{^-^)\l + 0{n-^)) ' ^ ' 

uniformly in arg A G [— f , ^i- 

Proof. Let A = -i(^)3(l + 0{n-^)) as n ^ +00. Then z = i^ + 0{n-^) or z = 
e~*f _l_ Consider the first case, the proof for the second one is similar. Using the 

identities 

e-i'^''2ei2(t. A) = e*(^-'^'')^e*(^"'^)^* = (^-(i+v^i)*)^^ 

' ' (7.4) 

e-'"'"ei3(t. A) = e'(^-"')"* = e""^"*, e-*"'"e23(t. A) = e'^"""')"* = e"^"* 

for all (t. A) G M X C, we obtain 

e-'"''ei2(t. A) = e-(^-')""*(l + ©(n^^)), e-*'"'"ei3(t. A) = e-(v^+^)""*(l + ©(n'^)). 



e"*'^ "e23(t. A) = e~*'™*(l + 0{n-^)) as n ^ +00 uniformly on t G [0, 1]. 

Substituting this asymptotics into identity ( 17. 2p and using ( 14. 9 p we obtain that the function 
ajk satisfies asymptotics (Q for A = -i{^f{l + 0{n-^)). 
Identities yield 

|e-*"''ei2(t, A)| = e-^("(^-*)+«*), |e-'"'^ei3(t, A)| = e~^^', |e-*"'^e23(t, A)| = e"^"* 

for all (t. A) G [0, 1] x C, where ^ = Substituting these estimates into identity (17. 2p we 

obtain 

|e-^'^'^a,fc(A)| ^ max [ f(u)g(u~t)du [ (e-^(^(^-*)+«*) + e"^^* + e-^^'*)rft (7.5) 
te[o,i] Jo Jo 

for all A G C. Let argA G [— f , Then argz G [— f^, ff] and max{x,^} ^ |2;|sinY^. 
Estimates (O]) yield (Q- ■ 

Lemma 7.2. Le^ p,q E L^(T). T/ien i/ie functions = Tr A^fc(l, ■), /c G N, satisfy 

T, = t^ T, = -^ + ^ + ^ T, = -i^ + ^ (76) 
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where $i is given by h2.21\) and 



e^{\) = J2 uj"'^e''^'\ m = 0,1,2, A G C, 

j=Q 



the functions T2,T3 satisfy ^7.3\ ), and 

Tfc(A) = e'^Od^r^) as \X\ oo, all k ^ 4. (7.7) 

Proof. Estimates fITTTD give ([721) • Identities ([221), O and C qit)dt = give 





1 2 







Then identities 



Ti(A) = Tr/ e^"(^-^)^Q(s,A)e'"*^rfs = Tre^^^ / Q(s, A)rfs = V e*"'^'"' / Qjj(s,A)rfs, 
Jo Jo Jo 

yield the first identity in (17. 6p . We will prove the second one. We have 



T2(A) = Tr/ / e''^^-'^^Q{t,X)e''^'-'^^Q{s,X)e'''^dsdt= y^^^jkW (7-9) 

where 

cijkiX) = / Cjkit - s,X)vjk{t,s,X)dsdt, Vjk{t, s, X) = Qjk{t, X)Qkj{s, X), (7.10) 
Jo ^/o 

and Cjfc has the form (17. ip . Identity (17.81) give 

a.,(A) = ^-(y^ Q,,(t,A)cit) =-^^^^^5 , all j = 1,2,3. 

Substituting these identities into (17. 9p we obtain 



T,{X) = -^^e,{X)+J2a,,{X). (7.11) 

Identity ([HD yields efcj(t, A) = ejk{l - t. A) for all j, k = 1,2, 3, {t. A) G M x C. Then fmU]) 
gives 

rl i*u i*\ 1*0 

o-kjiX) = du ejk{l — u + s, X)vjk{s,u, X)ds = dt ejk{t — u,X)vjk{t,u,X)du, 
Jo Jo Jo Jt-i 

which yields 

a.jfc(A) + afcj(A) = dt ejk{t - s, X)vjk(t, s, X)ds = / ejk{u. A) / vjkit, t - u, X)dtdu 
Jo Jt-i Jo Jo 

for all j, k = 1,2,3, X eC. Identities ([22D, give 
y^%fc(A)= V [\ue,k{u,X) [\,k{t,t-u,X)dt = ^^ + ^^ + 0(f^) (7.12) 
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as |A| — )■ cxD where 

$i(A) = -i [ (^jk{u,X) [ {p{t)q{t - u) + p{t - u)q{t))dtdu. 

By Lemma f7.1l the function $i satisfies f l7.3p . Substituting fl7.12p into (17. lip we obtain the 
second identity in (17. 6p and asymptotics (17. 3p for T2. 
We will prove identity (17. 6p for T3. We have 

1 t s 

T3(A) = Tr j j j e''^^-'^^Q{t,X)e''^'-'^^Q{s,X)e''^'-''^^Q{u,\)e''''^dudsdt 
000 

Identity (^M> and identity (EZD for Q give 

T3(A) = -^^^0 + ^( - 25i(A) + i?2(A)) + O(^) (7.13) 

as |A| — 7- 00 where 

3 3 
Bi = ^ uj^^'^^l3jk, I3jk = bjjk + bjkj + bkjj, B2 = ^ bj^^, (7.14) 
j,fc=i j,fc=i 

bjki= f f re^^("^''(^+"-*)+"'"(*-^)+-'"(^-"«/(t,s,«)^^«c^st^t, /(t,s,w)=p(t)p(s)p(w), 



Jo ./o Jo 

(7.15) 

7 = 7(j) ^) £ {1)2, 3}, 7 7^ j, 7 7^ A;. Assume that the functions -81,-82 satisfy asymptotics 
(I7.3p . Then identity (I7.13P gives the third identity in (17. 6p and asymptotics (17. 3p for T3. 

We will prove that the functions Bi, B2 satisfy asymptotics (17. 3p . Consider the function Bi. 
We have 

/•l l*t i*s i*\ nt 1*0 

^jjkW = / / / ejk{s — X)f{t, s,u)dudsdt = / / / ej^^t — X)f{t, s,u)dudsdt, 
Jo Jo Jo Jo Jo Jt-i 

fS 

bjkjW = / / / ejk{t — s, X)f{t, s,u)dudsdt, 
Jo Jo Jo 

1*1 i*t PS pi 1*0 i*s 

bkjjW = / / / ejk{l — t + u,X)f{t,s,u)dudsdt= / / / ejk{t — s,X)f{t,s,u)dudsdt. 
Jo Jo Jo Jo Jt-i Jt-i 

Substituting these identities into (I7.14p we obtain 

I t t 11 
Pjk{X) = j j j ejkit-s,X)f{t,s,u)dudsdt = j ejkiv,X) j p{t-v)p{t){p{t) -p{t-v))dtdv 

t-l s 

where p{t) = Jl^p{u)du. By Lemma [7.11 the functions I3jk, and then the function Bi, satisfy 
asymptotics (I7.3p . 
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Consider the function B2- Identities (I7.15P yield 



bjk-y{X) = dt ejk{t, s, X)h{t, s)ds, (7-16) 
Jo Jo 



bkj-yiX) = dt dsejk{t,s,X) j p{u — s)p{u — t)p{u)du 
Jo Jo 

for all A G C, j, 1 ^ j < A; ^ 3 where 



h{t,s)= / p{u - t + s)p{u - t)p{u)du, ejk{t,s,X) = ejk{t,X) 



Cjk are given by (17. ip . Consider the function 6123, the estimates for the other functions 
bjk-yyj, k = 1,2, 3,j 7^ k, are similar. We have 

e~^^^^ei2{t, s, X) = e«('^-'^^)2e*^(*-«-'^*+'^^^) = ^-^z{2-t-s-iV3it-s)) ^ (7.17) 

which gives |e-*^'^ei2(t, s,X) \ ^ e'^^^^"*) for all (t, s, X) e [0, 1]^ xC,s ^t. Substituting this 
estimate into (I7.16P and using the continuity of the function h{t, s) in s we obtain 

|e-*^'"6i23(A)| ^ [ dt [ e-^^(^-*)|/i(t,s)|c/s ^ [ max\h(t,s)\ [ e-^^'^^-'Usdt = O ( ^) 
Jo Jo Jo ^6[o.*l Jo ^FK 

as |A| — 00 uniformly in argA G [— f,^]- Thus the function 6123 satisfy (17.31) for argA G 

r 7r Stt I 

[ 4' 4 J- 

Let A = -i{^f{l+0{n-^)) as n ^ +00. Then z = i^+0{n-^) ot z = e"*! ^+0(n-^). 
Consider the first case, the proof for the second one is similar. Identities (I7.17P give 

e-'^^'^euit, s. A) = g-^^^^*"^) (e*™(*+^) + 0{n-^)). 

Substituting this asymptotics into (I7.16P and using (14. 9 p we obtain 



ie''"Hi23(A)| ^ dt e-^""(*-")|/i(t,s)|c/s(l + 0(n~i)) 
'0 Jo 



1 pI 

rfwe-^™" / \h{t,t-u)\dt{l + 0{n-^)) =0{n-^), 

J u 

which shows that the functions 6123 satisfy asymptotics (17. 3p . 

The similar arguments show that all functions bjk-y, and then the function B2, satisfy asymp- 
totics (17. 3p . which proves Lemma. ■ 

Proof of Lemma [Q Substituting (0, dUD, (EZD into the identity T = En>o^n 
obtain (^M>- ■ 
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